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Abstract 

A model is presented for generation of fast solar wind in coronal holes, relying on heating that is dominated 
by turbulent dissipation of MHD fluctuations transported upwards in the solar atmosphere. Scale-separated 
transport equations include large-scale fields, transverse Alfvenic fluctuations, and a small compressive dissi- 
pation due to parallel shears near the transition region. The model accounts for proton temperature, density, 
wind speed, and fluctuation amplitude as observed in remote sensing and in situ satellite data. 
Subject headings: MHD — waves — turbulence — solar wind 



1. INTRODUCTION 

An open question in solar and heliospheric physics is to 
identify the physical processes responsible for heating the 
corona and accelerating the fast solar wind streams emanating 
from coronal holes. This requires that a fraction of the energy 
available in photospheric motions be transported through the 
chromospheric transition region, and dissipated in the corona. 
The measured speeds of fast solar wind streams require spa- 
tially extended heating (Withbroe & Noyes 1977; Holzer & 
Leer 1980; Withbroe 1988). The physical mechanisms for 
this transport and dissipation have remained elusive. Some 
models have resorted to use of a parametrically defined heat 
deposition (a "heat function") that decays exponentially with 
height, or anomalous heat conduction that redistributes en- 
ergy along field-lines (Habbal et al. 1995; McKenzie et al. 
1995; Banaszkiewicz et al. 1998). One-dimensional (ID) 
models of this kind, extending from the chromosphere to 1 AU 
(Hansteen et al. 1994, 1997), have helped in understanding the 
regulation of the solar wind mass flux and can reproduce fast 
solar wind streams originating in cool electron coronal holes. 

Here we present a model that demonstrates solar wind ac- 
celeration due to heating by a quasi-incompressible turbulent 
cascade triggered by coronal stratification (Matthaeus et al. 
1999), and supplemented by compressive heating near the 
base. This model accounts for most presently available coro- 
nal and interplanetary observations. 

The idea that broadband plasma fluctuations might heat 
the extended corona and accelerate the solar wind has long 
been discussed (Coleman 1968; Belcher & Davis Jr. 1971; 
Hollweg 1986; Hollweg & Johnson 1988; Vefli 1993). How- 
ever, the mechanisms of transfer of fluctuation energy to small 
scales — and, in particular, the role of Alfvenic turbulence (as 
observed in the solar wind) and cascade processes (e.g., phase 
mixing, ponderomotive driving, shocks, etc.) — have not been 
described self-consistently. Two recent papers shed light on 
these relationships (Suzuki & Inutsuka 2005; Cranmer et al. 
2007). 
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Suzuki & Inutsuka (2005) incorporate ID compressive non- 
linear interactions driven by Alfven waves and leading to 
shock heating. This model produces good agreement with so- 
lar wind speed profiles. As low-frequency Alfven waves prop- 
agate upward, their wave pressure compresses the plasma. 
Unable to refract or mode-mode couple into a perpendicular 
wavenumber cascade, these waves must dissipate in ID shock 
fronts. This model provides a valuable demonstration that 
MHD fluctuations can act as a conduit to transport energy to 
the requisite altitudes. However the restriction to ID cascade 
is at odds with the well-established propensity for an incom- 
pressible MHD cascade to proceed mainly through wavevec- 
tors perpendicular to a strong mean magnetic field (Robinson 
& Rusbridge 1971; ShebaUn et al. 1983; Oughton et al. 1994; 
Bieber et al. 1996; Cho et al. 2002). Furthermore the corona 
exhibits a clear transverse structuring, and the initial fluctu- 
ations must have perpendicular correlation lengths not much 
larger than a super-granulation scale (15 000 km). 

Another recent model (Cranmer & van Ballegooijen 2005; 
Cranmer et al. 2007) incorporates a low-frequency cascade 
model (Verdini et al. 2006); however, the treatment of prop- 
agation and dissipation differs significantly from the present 
approach. Their scheme treats nonlinear effects as a pertur- 
bation, and it is unclear if it converges for strong turbulent 
heating. Here we employ a strong turbulence closure. We 
also do not rely on electron heat conduction to boost radial 
energy transport. Instead we compute an internal energy as- 
sociated with the protons only. This approach supports com- 
parisons with results employing improved representation of 
turbulence, such as shell models (Verdini et al. 2009a,b) and 
(potentially) full MHD simulation. 

We find here that reflection of Alfvenic turbulence alone 
does not lead to a full corona/solar wind stationary state - a 
compressible contribution is required. This is supported by 
recent observations from Hinode (Langangen et al. 2008; De- 
Pontieu et al. 2009). When we include a small component 
of compressive heating near the coronal base, fast solar wind 
streams are then accounted for. 



2. MODEL DESCRIPTION AND EQUATIONS 

We employ a ID large-scale steady-state MHD model in an 
expanding flux tube (cf. Leer et al. 1982; Verdini et al. 2006), 
comprising equations of mass continuity, radial momentum 
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conservation, and pressure (internal energy), 



to large-scale gradients by 



-[pf/A]=0, 



1 



U^^-ypV-U + {y-l)Q{r). 
ar 
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Here r is the radial coordinate, A(r) the flux tube cross- 
sectional area, p{r) - Ink-^T the thermal pressure, U - U{r)r 
the large-scale radial flow (wind) velocity, G the gravitational 
constant, and Msun the solar mass. R,- is the radial compo- 
nent of the (vector) divergence of the MHD Reynolds stress 
R = {Sb ■ V 6b I An — pSu ■ V6u}, where 6u, 6b are the fluc- 
tuations and the full magnetic field is B = rBr{r) + 6b, and 
f ■ 6b - Q - r ■ 6u. 
The total (thermal plus magnetic) pressure is 



6b^- 

P' ^Ink^T + . 



(4) 



We specify an area expansion factor A{r). Then Br{r) is de- 
termined by magnetic flux conservation, B,.(r)A(r) - const, 
which is constrained by 1 AU observations. 

Q{r) is the heating per unit volume, related to the heat 
function (per unit mass) H - Qjp. It involves an incom- 
pressible part Hi, associated with turbulence, modeled here in 
Karman-Taylor fashion (e.g., de Karman & Howarth 1938; 
Matthaeus et al. 2004). There is also a small compressive part 
He - Qclp, so that H - Hi + He- Turbulence influences the 
large-scale flow through Qi, wave pressure, and the Reynolds 
stress Rr, and is modeled using only a few free parameters. 

The dominant contributor to the turbulent heating, the low- 
frequency quasi-incompressible turbulence, is evolved us- 
ing a transport equation and one-point closure that depends 
upon the cross helicity. The nonlinear phenomenological 
model (Dmitruk et al. 2001) is involves the Elsasser variables 
- 6u + 6bl ^jAnp, their associated energies = Z±/4 - 
(|z^|^)/4, ((■■■) indicates an ensemble average) and a com- 
mon similarity (correlation) scale A. The dimensionless cross 
helicity cTc - (Zj -Z^)/{Zl +Z^) measures any excess inward 
or outward propagating-type fluctuations. The incompressible 
turbulent heating model(Hossain et al. 1995) is 



Hi(r) 



Qjjr) ^ l Z-Zl+Z+Z^_ 
p ~ 2 A 



(5) 



See also Dobrowolny et al. (1980), Grappin et al. (1983), and 
Matthaeus et al. (2004). 

We include spatial transport (Verdini & Velli 2007) in a 
non-uniform wind with speed U(r)r. The fluctuations are as- 
sumed to be Reduced MHD-like (i.e., perpendicular fluctu- 
ations, parallel gradients much weaker than transverse ones) 
and represented by "typical amplitudes" of a given frequency 



z±{(jS), defined such that Z| - ^ (z^/w^dw: 
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dr ^ 2A 
The WKB [Eq. (7)] and reflection [Eq. (8)] terms are related 
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Here Va = Va^ - B,-/ ^[^np. For simplicity the lengthscale A is 
associated with the expansion, using A(r) - A^ yfMr). The O 
subscript indicates evaluation at the coronal base, here taken 
to be at the top of the transition region. 

The smaller, compressive contribution to the heating {He) is 
assumed to be confined near the coronal base where it rapidly 
dissipates through shocks. This is motivated by recent obser- 
vations (Langangenet al. 2008; DePontieu et al. 2009) of fluc- 
tuations with parallel (vertical) variance that pervade the en- 
tire corona near the transition region. Here we model this ef- 
fect directly as a heat function, and assume that it contributes 
a! 1 % of the total (height-integrated) heat function. 

Having in hand the complete set of equations (l)-(3), and 
(6), along with constitutive relations (4), (5), (7), and (8), 
and the small compressive heat function He, we are in po- 
sition to solve a relatively complete and self-consistent solar 
wind model with turbulence. The stationary solutions for the 
large-scale fields and turbulence parameters are obtained via 
a numerical iteration procedure in which Tq = 4 x 10^ K and 
«0 = 5 X 10^ cm"^ are held constant. To begin, we impose 
trial radial profiles for T and n, and then solve Eqs. (6) with 
a fixed amplitude of the velocity fluctuation at the base, say 
6uq - 30kms"'. From this temporary solution we compute 
the Reynolds stress, the ponderomotive force, and the heating 
function that appear in Eqs. (l)-(3). Next, those equations are 
solved for T and U, which provides updated values for Va, 
U, and their derivatives. These updated profiles are used in 
re-solving Eqs. (6), and so on. 

After few iterations, there are small relative differences 
(< 10"^) in each of the large-scale and fluctuation fields. We 
control convergence through the relative variation of the posi- 
tions of the sonic and Alfvenic critical points, assuming con- 
vergence when Arj/rj, Arg/ra < 10""*. For each frequency 
a>, we solve subject to imposed conditions on 6uq, Hdr), the 
maximal (over)expansion /^ax - max[A(r)/r^], and the cor- 
relation length scale Aq. We use, as in Kopp & Holzer (1976), 
/(r) = A{r)/r, with 



f(r) 



/max exp[(r - r/)/g-] + fi 
exp[(r - r/)/cr] + /, 



(9) 



/i = 1 - /maxexp[(/?0 - r/)/cr], keeping fixed the location 
and width of the super-radial expansion, r/ - \3\Rq and 
cr - 0.5 1 Rq respectively. 

3. RESULTS 

A typical solution is shown in Figure 1 (solid black line), 
where we plot several quantities {U, T, n, (Je, Z~) as a 
function of distance, obtained with the chosen parameters; 
6uq = 30kms-i, w = OHz, ^ = 0.015 /^q, /^ax = 12.5, 
and He - \ %H. The speed U, temperature T, and density n 
are generally in good agreement with the observations in the 
heliosphere, although the terminal speed is about lOOkms ' 
higher and T peaks too far from the coronal base compared 
to the observed proton temperature. The latter is obtained 
by subtracting from the observed line width Wj. the contribu- 
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Fig. 1 . — Wind speed, temperature, density, turbulence level, and normalized cross helicity as a function of distance for the representative solution (solid line) 
with: (5i(o = 30 km s"', u) = OHz, Aq = 0.015 Rq^ /max = 12.5, and He = 1%H. The other lines denote solutions for which one of the parameters has been varied; 



specifically, 6uq = 20kms 



:0.05Rs,/„, 



10, and He ~ 5%H, in blue dotted, red dashed, green dashed-dotted, violet long-dashed lines respectively. 




Empty squares in the top-right panel mark the position of r, and (r, < r„). The other symbols represent observational constraints for the fast solar wind taken 
from: McComas et al. (2000) for U,n al I AU; Grail et al. (1996) for U inside 1 AU; Breech et al. (2005) for iTc and proton temperature beyond 1 AU; Cranmer 
(2009) for proton temperature inside 1 AU; Banerjee et al. (1998) and Fisher & Guhathakurta (1995) for n inside 1 AU; Bavassano et al. (2000) for the turbulence 
level. 

tion. This affects the turbulent energy Z^, which becomes 
larger beyond r„ x 13 Rq, and also 0-^., which stays closer to 
unity. U, T, and p are almost unchanged {U and T are slightly 
reduced). Below, we consider only w = fluctuations. 

The position of the sonic critical point depends upon mo- 
mentum and heat addition (Leer et al. 1982). Given that we 
specify Aq and p©, the mass flux is determined by Uq, which 
is found by requiring that the solution becomes supersonic on 
passing through r^. It follows that deposition of heat before 
increases T(rs) and typically Uq, and thus also the mass 
flux. Deposition of heat beyond r j does not alter T(rs), and 
so the mass flux is unchanged; however U increases and p de- 
creases beyond r,, with respect to a reference solution. Gen- 
erally speaking, the maximum of is below (it coincides 
with ra in the undamped case), the turbulent heating peaks 
near r^, and the momentum added by the wave pressure has a 
maximum inside r^. 

Two ways of controlling the heat and momentum deposition 
[i.e., adjusting the importance of nonlinear terms relative to 
reflection/WKB (linear) terms in Eqs. (6)], are through varia- 
tions of 6uq and A{r). Decreasing the basal fluctuation ampli- 
tude {Suq = 20kms"', in Fig. 1) decreases the turbulent heat- 
ing and acceleration. In addition heat deposition peaks well 
below rs (Fig. 2b), yielding a slow, overdense, cool wind, with 
a deficit of turbulent energy. It is clear that shapes the solu- 
tion only below 1 .5 Ro while differences arise from a different 
shape of //,. The asymptotic Z^ and cr^ do not vary much 
since beyond r,, reflection is negligible and the ratio Z_/Z+ is 
then controlled by turbulent dissipation and A{r). 

Modification of A{r) leaves unchanged the position of the 
critical points; i.e, it only slightly alters the turbulent heat- 
ing close to the coronal base, where reflection controls the 
ratio Z_/Z+. A larger basal correlation scale {Aq = 0.05 Rq in 
Fig. 1) affects the solution mainly beyond r,, yielding a hotter 
and faster wind, but with the density almost unchanged. The 
fluctuations have features similar to the undamped solutions, 
resulting in excess turbulent energy (that peaks closer to rg) 
and more inward propagating waves (a smaller cr^). Nonethe- 
less, the increased turbulent energy — driven by reflection and 
the WKB term — broadens the turbulent heating while keeping 
fixed its maximum (Fig. 2b). 

Adjusting the expansion A(r) also changes the correlation 
scale through A(r) 



Fig. 2. — (a). Total (thick line), compressive (thin line) and incompressible 
(dashed fine) heating function (dissipated energy per unit mass) as a func- 
tion of heliocentric distance in units of a 3 X lO'" cm^ s"' for the reference 
solution, (b). Total heating function for the solutions obtained varying the 
parameters as described in the caption of Fig. 1 (same line and color coding). 
Empty squares mark the location of r, for each solution. 



tion of the turbulent fluctuation 6u, according to the relation^ 



IkBT/m,, 



- {6uY/2 (e.g., Tu et al. 1998). The profiles of 



the fluctuation energy + - Z^/2 and normalized cross 
helicity ctc trend nicely towards the in situ data beyond 1 AU. 
The turbulent dissipation //, (Fig. 2a, dashed line) accounts 
for the spatially extended heating that accelerates the wind 
from below the sonic critical point. The compressive heating 
is small, about 1% of the total heating per unit mass. Its form 
is arbitrarily chosen to be a gaussian (Fig. 2a, thin solid line), 
centered at 1.3 Rq with width x 0.25 in order to confine 
its contribution well below rj. 

The solutions change as parameters are varied, but not 
all the parameters have the same impact. For example, for 
a steep spectrum of low-frequency waves (w < 10"^ Hz, 
slope < -1.1), use of only zero-frequency fluctuations is 
a very good approximation. For flatter spectra, the high- 
frequency part of the spectrum remains principally outward 
propagating [z^(aj)], thus limiting the total turbulent dissipa- 

' Only the reference case is used to compute the corrected proton temper- 
ature in Fig. 1. Other cases give similar values except the 6uq = 20kms"' 
and He = 5%H cases, yielding much higher coiTected temperatures. 



Aq yA{r). A(r) controls the reflec- 
tion term directly and through the Alfven speed gradients 
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[Eqns. (7)-(8)]. Decreasing the maximum super-radial ex- 
pansion, to /max = 10, while keeping fixed its location at 
Vf Si 1.3^0 < (Fig. 1) has two consequences. First, it 
decreases the correlation scale for r > rj, causing a reduc- 
tion of turbulent energy and heating. Second, it increases the 
density scale height around r/, hence reducing reflection and 
the amount of turbulent dissipation there. The result is much 
smaller turbulent heating that peaks again at (Fig. 2b) yield- 
ing a slower, cooler, and denser wind, even as close as r^, 
which retains good asymptotic and cr^. 

Finally, let us examine the role of compressive heating, re- 
calling that H{r) - He + //,. If the compressive heating is in- 
creased to He ~ 5%H (Fig. 1, violet dashed), the temperature 
maximum moves closer to the coronal base, but the resulting 
wind is slower, denser, and cooler An increased H^ also alters 
the incompressible heating — the height-integrated heating is 
almost unchanged, but heat deposition occurs at lower r due 
to increased reflection at r < rj (Fig. 2b, compare the smaller 
(Tc in Fig. 1). Beyond the peak of H the fluctuation energy 
remains small, reducing the extended turbulent heating. 

On the other hand, decreasing He enhances the incompress- 
ible heating; then H has a minimum near the maximum of 
Va, just inside r,, and also peaks outside r,. This causes 
convergence problems. In early iterations, heat deposition is 
mainly outside r,, yielding a very fast, underdense, and hot 
wind, with strong density gradient (and hence reflection) at 
r > Ts- This leads to runaway iterations in which the sonic 
point moves to larger r with a temperature minimum inside 
rs- When the temperature minimum becomes very low the 
sonic critical point is reached with dT/dr < 0. Then steep 
gradients occur close to the coronal base, heat deposition oc- 
curs close to r,, and the solution again resembles the initial 
solution, producing the ensuing runaway. To obtain a solu- 
tion, the key features of He are that it is localized well below 
r,, and is at least ~ 0.7% of the total heating (for other param- 
eters at reasonable levels). We found that the results do not 
depend strongly on the form of the compressive heat function 
provided that this heating is localized and not too large (expo- 
nential, gaussian, etc. all work). 

4. DISCUSSION 

The above model shows that turbulence near the coronal 
base, originating through chromospheric transmission of fluc- 
tuations, can heat the plasma in an expanding coronal hole 
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flux-tube and produce a fast solar wind that matches a num- 
ber of observational constraints. The turbulence is mainly of 
the low-frequency Alfvenic type. A small amount of com- 
pressive heating between the transition region and the sonic 
point appears to be needed to match the observations. This 
additional heating may be due to type II spicules that supply 
broadband low-frequency vertical fluctuations at transition re- 
gion heights, thus launching compressive MHD modes near 
the coronal base (DePontieu et al. 2009). 

Most of the fluctuation energy is in low-frequency turbu- 
lence, and this sustains a strong anisotropic MHD cascade 
through reflections from local density gradients. This type 
of anisotropic cascade is favored in MHD turbulence in the 
presence of a strong DC magnetic field (Robinson & Rus- 
bridge 1971; ShebaHn et al. 1983; Oughton et al. 1994). Heat 
conduction does not enter the present model at all, since it 
mainly affects electron internal energy, which evolves inde- 
pendently in this approximation. Similar assumptions work 
well in understanding observations of solar wind turbulence 
(Breech et al. 2009; Cranmer et al. 2009). In these ways, 
the present model differs from other recent models that incor- 
porate turbulent heating (Suzuki & Inutsuka 2005; Cranmer 
et al. 2007). In particular we believe that this model demon- 
strates, possibly for the first time, that a model (almost) free 
of ad hoc heat functions, artificial equations of state, and ad 
hoc assumptions about heat conduction, can indeed heat the 
corona and accelerate the solar wind. 

We plan further study of the required small amount of com- 
pressible heating, attributed here to spicule-driven magne- 
tosonic activity. Another useful extension would be to include 
separate electron and proton internal energy budgets, which 
will enable additional observational constraints, and will per- 
mit study of the role of kinetic dissipation processes (Breech 
et al. 2009; Cranmer et al. 2009). 



This research was supported in part by NASA (Helio- 
physics Theory, NNX08AI47G) and NSF (Solar-Terrestrial 
and SHINE programs ATM-0539995 and ATM-0752135). 
AV acknowledges support from the Belgian Federal Science 
Policy Office through the ESA-PRODEX program. We thank 
the UVCS project (J. Kohl, S. Cranmer) for providing data 
used in Fig. 1. 



DePontieu, B., Mcintosh, S. W., Hansteen, V. H., & Schrijver, C. J. 2009, 
ApJ, 701, LI 

Dmitruk, P, Milano, L. J., & Matthaeus, W. H. 2001, ApJ, 548, 482 

Dobrowolny, M., Mangeney, A., & Veltri, P. 1980, Phys. Rev. Lett., 45, 144 

Fisher, R. & Guhathakurta, M. 1995, ApJ, 447, L139+ 

Grall, R. R., Coles, W. A., Klinglesmith, M. T., et al. 1996, Nature, 379, 429 

Grappin, R., Leorat, J., & Pouquet, A. 1983, A&A, 126, 51 

Habbal, S. R., Esser, R., Guhathakurta, M., & Fisher, R. R. 1995, 

Geophys. Res. Lett., 22, 1465 
Hansteen, V. H., Leer, E., & Holzer, T. E. 1994, ApJ, 428, 843 
Hansteen, V. H., Leer, E., & Holzer, T. E. 1997, ApJ, 482, 498 
Hollweg, J. V. 1986, J. Geophys. Res., 91, 4111 
Hollweg, J. V. & Johnson, W. 1988, J. Geophys. Res., 93, 9547 
Holzer, T. & Leer, E. 1980, J. Geophys. Res., 85, 4665 
Hossain, M., Gray, P. C, Pontius, D. H., Matthaeus, W. H., & Oughton, S. 

1995, Phys. Fluids, 7, 2886 
Kopp, R. A. & Holzer, T. E. 1976, Solar Phys., 49, 43 
Langangen, 0., De Pontieu, B., Carlsson, M., et al. 2008, ApJ, 679, L167 
Leer, E., Holzer, T. E., & Fla, T. 1982, Space Science Reviews, 33, 161 
Matthaeus, W. H., Minnie, J., Breech, B., et al. 2004, Geophys. Res. Lett., 

31, 12803 



Turbulence-driven fast solar wind 



5 



Matthaeus, W. H., Zank, G. P., Oughton, S., Mullan, D. J., & Dmitruk, P. 

1999, ApJ. 523, L93 
McComas, D. J. et al. 2000, J. Geophys. Res., 105, 10419 
McKenzie, J. R, Banaszkiewicz, M., & Axford, W. I. 1995, A&A, 303, L45 
Oughton, S., Priest, E. R., & Matthaeus, W. H. 1994, J. Fluid Mech., 280, 95 
Robinson, D. C. & Rusbridge, M. G. 1971, Phys. Fluids, 14, 2499 
Shebalin, J. V., Matthaeus, W. H., & Montgomery, D. 1983, J. Plasma Phys., 

29, 525 

Suzuki, T. K. & Inutsuka, S. 2005, ApJ, 632, L49 

Tu, C.-Y, Marsch, E., Wilhelm, K., & Curdt, W. 1998, ApJ, 503, 475 

VeUi, M. 1993, Adv. Space Res., 14, 123 



Verdini, A., Dmitruk, P, Matthaeus, W., Oughton, S., & Velli, M. 2006, 

SOHO-17. 10 Years of SOHO and Beyond, 617, 150 
Verdini, A. & Velli, M. 2007, ApJ, 662, 669 

Verdini, A., Velli, M., & Buchlin, E. 2009a, Earth Moon and Planets, 104, 
121 

Verdini, A., Velli, M., & Buchlin, E. 2009b, ApJ, 700, L39 

Withbroe, G. L. 1988, ApJ, 325, 442 

Withbroe, G. L. & Noyes, R. W. 1977, ARA&A, 15, 363 



